The physical picture behind the oscillating sign of drag in high Landau levels 
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We consider the oscillating sign of the drag resistivity and its anomalous temperature dependence 
discovered experimentally in a bi-layer system in the regime of the integer quantum Hall effect. We 
attribute the oscillating sign to the effect of disorder on the relation between an adiabatic momentum 
transfer to an electron and the displacement of its position. While in the absence of any Landau 
level mixing a momentum transfer ftq implies a displacement of ql 2 H (with Ih being the magnetic 
length), Landau level mixing induced by short range disorder adds a potentially large displacement 
that depends on the electron's energy, with the sign being odd with respect to the distance of that 
energy from the center of the Landau level. We show how the oscillating sign of drag disappears 
when the disorder is smooth and when the electronic states are localized. 
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Coulomb drag in bi-layer systems at high Landau levels 
has been an active area of research in the last few years, 
both theoretically and experimentally. This research has 
been initiated by the puzzling results obtained in a series 
of experiments 2 **^, where the longitudinal drag resistiv- 
ity p D was measured in the presence of strong magnetic 
fields and large Landau level filling factors. 

In a typical drag experiment, a current j( a ' is driven 
through one layer, the active layer, while no current is 
permitted to flow in the other layer, the passive layer. 
The longitudinal drag resistivity is defined by 
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where is the component of the induced electric field 
in the passive layer that is parallel to p a \ 

Experimentally, two anomalous features of p D at low 
temperatures were discovered. The first is a sign that 
oscillates when the density difference between the two 
layers is varied: when the Landau level filling factors of 
the two layers differ by an even number p D is positive, 
i.e., has the same sign as at zero magnetic field. In con- 
trast, it becomes negative when the filling factors of the 
two layers differ by an odd number. The second anomaly 
is the non monotonous dependence of p D on the tem- 
perature T. For high temperatures, p D has a quadratic 
dependence on T, similar to that observed at zero mag- 
netic field. However, as temperature is decreased to low 
enough values there is a sudden increase of | p" 1 . As the 
temperature is further decreased \p D \ decreases, finally 
reaching zero in the limit of T — * 0. Both anomalous 
features occur at the same range of temperatures. 

The experimental observations were theoretically ad- 
dressed in several works^i^. An important step to- 
wards their understanding has been carried out by Gornyi 
et al£, who used the Self Consistent Born Approxima- 
tion (SCBA), and found the sign of p D to oscillate as the 
densities of the layers are varied. The non-monotonous 
temperature dependence of p D was obtained as well, al- 
beit its Arrhenius functional form at low temperature 
was missed. Although the SCBA calculation is generally 
consistent with the experimental observation, the physi- 
cal picture behind its results is obscure. 



In this work we study p D in a way that makes this 
physical picture transparent. We show that the oscillat- 
ing sign of p D originates from the effect of disorder on the 
relation between an adiabatic momentum transfer to an 
electron and the displacement of its position. Coulomb 
drag stems from scattering processes in which an energy 
lu < T and a momentum q are transferred from one layer 
to another (we set h = 1 throughout this paper). In 
the limit of low temperature, u> is very small, and the 
transfer of momentum is adiabatic. When disorder is 
weak enough to allow any Landau level mixing to be ne- 
glected, this momentum transfer implies a displacement 
of ql 2 H . We find that with the inclusion of lowest order 
Landau level mixing, the combined effect of a short range 
disorder and the rapidly oscillating wave functions char- 
acteristic of high Landau levels results in an additional 
displacement. The sign of this displacement depends on 
the electron's energy relative to the center of its Landau 
level, and its magnitude becomes large in the limit of 
a high Landau level. The electron's energy is typically 
the chemical potential, and thus p D develops oscillations 
with respect to the variation of the chemical potential in 
either one of the layers. Furthermore, we show that the 
non-monotonic temperature dependence originates from 
the strong oscillations in the density of states character- 
istic of well separated Landau levels. 

We consider the weak coupling regimei2*iii^ in which 
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where py$ and piy are respectively the Hall resistivities 
of the passive and active layers, U is the screened inter- 
layer interaction and B is the magnetic field. Most im- 
portantly, r is the rectification function relating a scalar 
potential Vq U of wave vector q and frequency to to the 
DC current it creates in a second order response: 
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The function T is a vector, with a longitudinal com- 
ponent r 1 1 parallel to the wave vector q and a transverse 
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component T± perpendicular to q. For symmetry rea- 
sons, r satisfies 

r||(q,w,B) - _r|,(-q,«,-B) (4) 
rj.(q,w,B) - rj.(-q,w,-B) (5) 
r(-q,w,B) = -r(q,w,B). (6) 

Eq. © is valid provided T is averaged over disorder. 

Given Eqs.|[2HfjJ and the experimental observation of 
positive p D for identical layers, it is clear that the drag 
resistivity is dominated by the contribution of r_i_. Thus, 
the oscillations of the sign of p D and its T dependence 
should originate from the dependence of T± on u, q and 
the Landau level filling factor (or, equivalently, the chem- 
ical potential). In particular, the density dependence of 
the sign of p D must originate from a dependence of the 
sign of T± on the chemical potential^. 

In view of Eq. © we find it more convenient to define 
the Hall rectification function 

r H (*h u) = \ (r±(q, uj, B) + r ± (- q , «, —B) \ . (7) 

We consider a model of non-interacting electrons sub- 
jected to a magnetic field B and a disorder potential Vd- 
The single layer Hamiltonian is then, 

tf = £{^(P.- A M) 2 + ^r;)}. (8) 
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This Hamiltonian can, in principle, be diagonalized by 
single particle states with energies {e Q } and wave func- 
tions {"fa}- Hence, the unperturbed Green functions ob- 
tain a simple single particle form, and in the basis of 
exact single particle states Th is given by the triangle 
diagramsiiiii. 

As we show below, within this model, 

r H (q,u>) = Y / Ax0 J T 0r . (9) 

0-y 

where 

Tf3 7 = 27r|(p q ) 7/3 | 2 [n F (e 7 ) - n F (ep)] 8(ep - oj - e 7 ) (10) 
and 

Ax fa = {x)f3/3 - (a;) 77 . (11) 

Hereafter the y direction is set parallel to q, p q is the 
density operator, n F is the Fermi distribution function, 
and {x)f3f3 is the expectation value of the position x with 
respect to the single particle state We discuss below 
the way Eq.® is to be understood for extended states, 
for which this expectation value may be ill defined. 

Eq.© gives a very simple physical explanation to the 
Hall rectification current r#(q, w)|Vq W | 2 . The applica- 
tion of a time and space dependent potential Vq W induces 
transitions between single particle states. The net rate 
for transitions between the states 7 and f3 is given by 



Fermi's golden rule to be T , ( 3 7 |Vq W | 2 . A transition from 
a state 7 to a state j3 involves also a translation of the 
electron's position given by Axp^ , and thus induces a 
current. The total current involves a sum over all possi- 
ble transitions, and hence Eq.©. 

Furthermore, Eq.© clarifies that the sign of Th is de- 
termined by the preferred direction of the induced transi- 
tions. Put differently, to determine that sign one should 
identify the states (3 and 7 that make the most signifi- 
cant contribution to TJg 7 , and determine the sign of the 
displacement Axpy associated with their transition. 

With the imaginary part of the polarization operator 
being Imn = X)/3 7 ^87 we nn( i it instructive below to 
express r h in terms of Imn. One can view the ratio 
between them as the characteristic displacement for a 
particular q, uj and chemical potential. 

In our discussion we are guided by the experimental 
conditions to focus on well separated Landau levels, large 
filling factors and q <C fc F , with k F being the Fermi mo- 
mentum. We denote the uppermost filled Landau level 
by n, and consider the limit of low temperatures. We also 
find it convenient to measure all energies with respect to 
(n + l/2)w c , the center of the n'th Landau level. 

We first discuss the expectation value {x)pp. The sin- 
gle particle states \/3) may be written as |/3o) + |/?i); where 
\/3o) is composed of states of the n'th Landau level, and 
|/3i) is composed of states of other Landau levels. To lead- 
ing order in the disorder potential, \/3q) — ^2kC^\nk), 
where {|nA:)} are the n'th Landau level wave functions of 
a clean system in the Landau gauge, A = B(0,x). The 
coefficients and the leading order approximation to 
the energies are found within first order degenerate 
perturbation theory to satisfy, 

EHI^M^ = (12) 

p 

The Landau level mixing part, is calculated by 

means of first order non-degenerate perturbation theory 
and is of first order in the ratio Vd/oj c . Within that order, 

(aO/s/s = (A)M#>) -iu}- 2 (p \[v x ,V d ]\f3o). (13) 
where we use the decomposition 

v 

r = r - z x — (14) 

of the position operator r, with ro the guiding center co- 
ordinate and v the velocity operator. In Eq. (|13|l the 
first term is the expectation value of the guiding cen- 
ter coordinate xq, and the second term is the expecta- 
tion value of Vy/io c , which becomes non-zero only due to 
the Landau level mixing induced by disorder. Note that 
[v x ,Vd(x)] oc d x Vd{x). The shift of the position induced 
by the disorder is proportional to the expectation value 
of the electric field experienced by the electron. 

We now turn to discuss the states \(3o), |7o)- Within 
the semiclassical approximation, 

(7oMA)> = MqRc)ho\e- iqv °\f3 ) (15) 
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with R c being the cyclotron radius of the n'th Landau 
level. Thus, in order to analyze the most significant 
transitions, we consider two states: the first is a sin- 
gle particle eigenstate |/3q) with an energy ep, in the 
vicinity of the chemical potential fi. The second state, 
70} = e~~ lqvo \Po), is a state defined in such a way that 
the potential V(q,u>) couples it most effectively to |/?o) : 
its momentum is shifted by q relative to that of |/3o)j an d 
it is projected to the n'th Landau level. Furthermore, 
focusing on the low temperature limit of Eqs. (|2I1U|) we 
look for those cases for which the expectation value of 
the energy of I70) equals tp, and find the implications of 
this condition on Aa: / 3 7 . 

When the chemical potential is in a region of local- 
ized states, |/3o) and I70) are both localized, and the ex- 
pectation value of the velocity, the second term in 1(13(1 . 
vanishes. Therefore, 

= (^Wo - (e^zoe-^Wo = ql 2 H , (16) 

and 

T% oc) = 2< 7 Z 2 f ImII qa) . (17) 

The sign of Th is then independent of p. 

When the states |/?o), |7o) are extended the second 
term in 1(13(1 no longer vanishes and 

^4f=ql 2 H + ^Mq)-v (O)), (18) 

where vp(q) = d q e((3,q) with 

e(0,q) = J dv\^ n p{x~ql 2 H ,y)\ 2 V d (x,y) (19) 

being the energy of the state e~ tqy ° \(3o). For a slowly 
varying potential changing on a scale of £ > R c the sec- 
ond term of Eq. I|18|) is of the order of ql 2 H which is 

negligible compared to ql 2 H . When that happens relation 
(|17fl between Th and Imll still holds, although Imll now 
reflects the extended nature of the single particle states. 
The opposite extreme is exemplified by 

V d (r)=^26(r-r } )U j (20) 
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describing point like impurities. We assume the impuri- 
ties' positions {r,,} to be random, and their strengths Uj 
to be symmetric around zero. For the simplicity of the 
calculations, we use a rather crude approximation for the 
\nk) states, 

piky 

\nk) = cos(k F (x-kl 2 H ))Q Rc (x-kl 2 H ) (21) 

V RcLy 

where L y is the system's length in the y direction and 
Qr(x) is a smoothed step function that equals one for 
x 2 < R 2 and zero otherwise. For any fixed x, Eq. I|21|) 



becomes exact in the limit n — * 00. Within the approxi- 
mation (|21|l . the Bessel fun ction in 1)15(1 is replaced by a 
cosine. It follows from Eas. l|20l21|l that 

e(0,q) = (RcLy)- 1 ^ C^Cfae-^-*'^ cos(e jp + qR c ) 
p P 'j 

cos(6» jy +qR c )Q Rc (x ] ~pl 2 H )Q Rc {x 3 -p'4),(22) 

where 9j p — k F (xj —pl^). Since q -C k F we neglected the 
change in the argument of the 9 functions . Making the 
assumption that for all y, the disorder potential averages 
to zero when integrated over x, i.e., J dxVd{x, y) — 0, we 
obtain from Ea. 1(22(1 a relation between the velocity vp(0) 
and the energies ep and e 7 : 

e 7 = cos(2qR c )ep - 7TFr v p(°) sin(2qR c ). (23) 
ZH C 

Being interested in the low frequency limit of 1(10(1 . we 
look for states |/3o) for which ep = t 1 . For these states, 

v p (0) = -2R c eptim(qR c ). (24) 

At low temperatures ep ~ \i and hence the sign of f/3(0) 
depends on the position of fi relative to the center of the 
Landau level. By interchanging the role of states (3q and 
70 we find vp(q) = —vp(Q). Thus, 

T% xt) a (2ql% + 4i? c £- tan(qR c ^j ImU^. (25) 

The second term in the brackets is the crucial disordered- 
induced contribution to the displacement, whose sign 
changes as n crosses the center of the Landau level. This 
contribution, which stems from the second term in 118II . 
becomes larger than the first term for a certain range of 
q. Although it is sub- leading in Vd/ui c , this range is made 
large at high Landau levels due to the dependence of 1(24(1 
on i? c , a dependence that originates from the rapid os- 
cillations of the wave functions on the scale of the Fermi 
wavelength. When this range is large enough, the change 
of sign of r# as a function of fi leads to the sign vari- 
ations of p D &. We note that within approximation 1(21(1 

Imn qw oc cos 2 (qR c ) so that Y is finite for all qR c . 

Ea. ((25)) is similar to the SCBA result^. The two re- 
sults differ only in the replacement of Bessel functions by 
trigonometric ones. The discrepancy may result from ap- 
proximation ((21(1 that neglected the increase in the wave 
function amplitude as the classical turning points are ap- 
proached. 

The extended states lie at the center of the Landau 
level whereas the localized states lie at its ends. Thus, if 
the chemical potential lies in the region of the localized 
states, the contribution of r^*-* will show activated be- 
havior at low temperatures. The relative weight of the 
contributions of localized and extended states is difficult 
to estimate, but it is possible that this activated behavior 
is the source of that seen in the experiments. 

The above physical understanding of Th sheds light 
on the non monotonous temperature dependence of p D . 
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As shown in Th is a weighted sum over all possible 
single electron transitions. The transition strength T ( g 7 , 
which is a measure of the net transition rate between the 
two states, is non zero only if ep = e 7 + uj. From Eq.J2J 
it follows that uj < T. For well separated Landau lev- 
els, A <C uj c , with A being the Landau level's width. At 
low temperatures, T <C A, the phase space for scattering 
increases with temperature. Hence, Th and p D increase 
with temperature as well. However, for A < T < w c 
the difference in occupations n F (e 7 ) — n F (ep) scales as 
1/T and decreases with increasing temperature. Con- 
sequently so does the net transition rate between these 
states, as well as Th and p D . Finally, for T ^ oj c inter 
Landau level transitions become possible and p D returns 
to increase with temperature. To fully account for the 
temperature dependence of p D the temperature depen- 
dence of the screened inter-layer interaction, U , must be 
incorporated as well£. Nevertheless, the above physical 
picture remains intact. 

We now turn to the derivation of Eq.@. In terms of 
the exact eigenstates {^ a } and exact energies {e a }, 

^{v x ) ia {p q ) a p{P- q )M 

xF(e a ,e /3 ,e 7 ,w), (26) 

where F is determined by the standard Keldysh 
techniqua2ii2iii. The inversion of the magnetic 
field leaves the energies unchanged while transforming 
{$>a} — > {^*a\- Thus, the second term of T H , given by 
Eq. , can be related to the first term. As evident from 
Eq.© the rectification function is real. Hence 

— Im{(u x ) 7a (/3 q ) Q/3 (p_ q ) / 3 7 } 

x Re { G r a (e) (G£ (e + w) - G% (e + uj)) G« (c) } 

x (n P (e)-n P (e + w)) + \ * L (27) 

where G r a and G° are, respectively, the retarded and 
advanced Green functions. After some straightfor- 
ward manipulations of which we particularly note the 



use of Heisenberg's equation of motion, (v x ) ia 
i (e 7 — e Q ) (x) 7 a, we obtain 

r H (q,u) = 27r^ (n F (e 7 ) - n F (ep)) 6(ep - co - e 7 ) 

0i 



Re { (3^)7/3 (P-q)/37 - (x) 77 |(p q ) 7/3 | 2 } + 




(28) 

Adding the q — » — q, w — * — lo term explicitly and using 
[n;,pq] = we obtain Eq.©. 

Ea . (|28fl is valid provided that As^ 7 is finite. This 
is certainly true for localized states. It follows from 
Ecis. l|14ll6jl that Ax^ 7 is finite for extended states as well 
since (2:0)77 — {xo)pf3 ~ Qlffi despite the fact that both 
of these expectation values are ill defined. 

In summary, we studied the drag resistivity, p D , for 
well separated high Landau levels, focusing on the de- 
pendence of the Hall rectification function T# on the 
chemical potential p, and temperature T. Using the ba- 
sis of exact single particle eigenstates, we expressed Th 
as a sum over transitions induced by the potential Vq W . 
The contribution of each transition is the product of its 
rate by the displacement it induces. For localized states 
the displacement is ql H , yielding a sign of Th that is in- 
dependent of p. For extended states this displacement 
is augmented by another contribution, induced by the 
disorder potential. The sign of the latter contribution 
is odd with respect to the distance of p from the cen- 
ter of the Landau level. For high enough Landau levels 
this contribution dominates and induces oscillations in 
Pd as a function of density difference between the layers. 
The non monotonous dependence of p D on T is a conse- 
quence of the strong oscillations in the density of states 
of well separated Landau levels. The difference between 
r# for localized and extended states may account for the 
activated T dependence of p D . Our results for Th of ex- 
tended states with short range disorder are in agreement 
with the SCBA results of Gornyi et al£, and unravel the 
physical picture behind this approximation. 
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